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We study the surface critical behavior of semi-infinite quenched random Ising-like systems at the 
special transition using three dimensional massive field theory up to the two-loop approximation. 
Besides, we extend up to the next-to leading order, the previous first-order results of the ex- 
pansion obtained by Ohno and Okabe [Phys. Rev. B 46, 5917 (1992)]. The numerical estimates 
for surface critical exponents in both cases are computed by means of the Fade analysis. Moreover, 
in the case of the massive field theory we perform Pade-Borel resummation of the resulting two- 
loop series expansions for surface critical exponents. The obtained results confirm that in a system 
with quenched bulk randomness the plane boundary is characterized by a new set of surface critical 
exponents. 



I. INTRODUCTION 



The investigation of the critical behavior of the real physical systems is of considerable theoretical and experimental 
interest. As usual, the real physical systems are characterized by the presence of different kinds of imperfections. In 
common, the defects and impurities may be localized inside the bulk as well as at the boundary. 

Historically the systematic investigation of the quenched disordered systems was initiated in the seminal works 
by Harris, Lubensky [^]|^] and Khmelnitskii It is shown that investigation of the Ising-like systems is of main 
concern from the whole class of 0{n) symmetric n-vector models in d dimensions, because they satisfy the Harris 
criterion Q|, which states that the presence of randomness is relevant for such pure systems which have a positive 
specific heat exponent a. The introducing of the bulk dilution into a system shifts the critical temperature of the 
bulk phase transition and drives the system to another, 'random' fixed point in which unconventional scaling behavior 
is observed. As it is confirmed by the Wilson's renormalization group and e expansions |[l[-^,^-Q , the massive field 
theory in three dimensions , experiments ||l^-|l^ , and Monte-Carlo simulations ||l^;|^ , the critical behavior of 

three-dimensional disordered Ising-like systems is characterized by a new set of critical exponents |l7|. The case of 
the Ising model at d = 2 is marginal, because in this case a — and the correspondent logarithmic corrections to the 
power laws singularities of the pure model take place, as was confirmed in a series of papers p0|-|25[|. 

The presence of a surface leads to the appearance of additional problems in critical phenomena. The most general 
review of critical behavior at surfaces and list of related publications are given in [p6|-p8[ . It is well known p^-pO|] that 
each surface universality class is defined by the bulk universality class and specific properties of a given boundary. 
At the present time three surface universality classes, called ordinary, special and extraordinary, are known p7| , ^ . 
They correspond to the respective surface transitions which occur at the bulk critical point r7i§ — tuq^ and are 
characterized by different fixed points 

*,ord , sp *,extr /1 1 \ 

Cq' = +00, Co' = Cq"^, Cq' = -oo. (1.1) 

Here cq is so called 'bare surface enhancement', which measures the enhancement of the interactions at the surface, 
and (mQ,co) = {mQ^,CQ^) is a multicritical point, called special point. 

The influence of quenched surface disorder on the surface critical behavior was investigated in a series of theoretical 
works 1^,^ and Monter Carlo calculations [^,^. General irrelevance-relevance criterion of the Harris type for 
short-range as well as for l ong-ra nge correlated surface disorder was derived In the case of special transition 

it has been demonstrated ||3^,^ that the fixed point describing the surface critical behavior of three dimensional 
pure systems is stable with respect to short-range correlated surface disorder. Thus, the weak short-range surface 
disorder is irrelevant for three dimensional systems, but long-range correlated random field disorder is relevant in 
d < 4 dimensions. 

In a recent paper js^, we quantitatively confirm the previous results by Ohno and Okabe p5|| that the introducing 
of quenched bulk randomness in semi-infinite systems bounded by a plane surface affects the surface critical behavior 
of these systems. Thus obtained surface critical exponents of quenched dilute semi-infinite systems at the ordinary 
transition differ from the surface critical exponents of the pure semi-infinite systems p5| . Besides, we showed that to 
order e, the ^/e expansion for surface critical exponents t^h and rij_ has given negative value of the correlation function 
critical exponent r/ for the random bulk Ising system according to the scaling relation rj = 2ri± — . It confirms 



1 



the well known fact that the second order of thcy^ expansion is not enough to give correct positive value of critical 
exponent rj P,|7|,|37|-^ . The obtained results |Q have shown that this kind of deficiencies do not appear at the 
calculations using the massive field-theoretic approach directly in d = 3 dimensions [^ . 

All these have stimulated us to perform the investigation of the special transition occurring in quenched bulk dilute 
semi-infinite systems bounded with a plane surface. It should be mentioned that the problem of investigation of the 
special transition is very important from such point of view that at some conditions it may be reduced to the problem 
of the adsorption of 9 polymers on a wall . 

The investigation of the critical behavior of the systems with quenched, i.e., time independent, randomness is 
possible to split on two directions. One of them is renormalization-group approach introduced by Harris and Lubensky 
[01 . This approach involves applying the renormalization-group transformation to the random system directly and 
subsequent averaging over disorder. Ohno and Okabe |4q| employed the above mentioned method to analyze the 
influence of randomness on the surface critical behavior at d = 4 — e dimensions in the frames of ^/e expansion. 

Another technique introduced by Grinstein and Luther |^] involves firstly removing the randomness by averaging, 
and subsequent employing the renormalization group. They considered an mn vector model and showed that analytic 
continuation of this model to rt = is equivalent to a model of a random m-component spin system. An elegant 
derivation of this equivalence has been given by Emery p6|| . Our main calculations are performed with this technique 
to treat randomness. 

The present paper is dedicated to the investigation of the special surface transition in semi-infinite, quenched dilute 
Ising-like systems at the bulk 'random' critical point directly in c? = 3 dimensions using the massive field theory up 
to the two-loop approximation. Besides, we extend up to the next-to leading order of the ^/e expansion, the previous 
first-order results obtained by Ohno and Okabe p5[ . The numerical estimates for surface critical exponents of the 
special transition in both cases are calculated using extensive Pade analyses. Moreover, in the case of the massive 
field theory we performed Pade - Borel resummations of the resulting two-loop series expansions and obtained quite 
reasonable and reliable numerical estimates for surface critical exponents. The obtained results confirm that in the 
case of quenched bulk randomness in semi- infinite systems the new set of the surface critical exponent appears. 



II. MODEL 

In the previous work p6[ |, we presented some notations about possibilities to use effective Landau-Ginzburg- Wilson 
Hamiltonian with cubic anisotropy defined in semi-infinite space 







Hlgw\<P] = I dzld'^-^r 



i|V^P + imgl'?l' 



Q=l 



ivo E <l>t + iuo{m' , (2.1) 



for description of surface critical behavior of quenched dilute semi-infinite Ising-like systems at the ordinary transition 
in the replica limit n — > 0. It should be mentioned that here (/) is an n - vector field with the components 4n,i = \, ...,n 

defined on a half-space M'\. = {x=(r, z) e M'^ \ r e M'^~^,z > 0} bounded by a plane free surface at 2; = 0. 

For the first time the above model was introduced by Grinstein and Luther |^|. The 0(n) symmetric term in (2.1) 
arises in the process of the configurational averaging of the free energy over disorder using the replica trick 

F = -Tlim-{{Z"),,^f-l), (2.2) 

n^O n 

via cumulant expansion and in accordance with this its coupling constant uq c< < 0- Here Z is the partition 
function with Boltzman weight e^^^"^! , where H[(j)] is the effective Landau-Ginzburg- Wilson Hamiltonian of an original 
Ising system with scalar field — 4>{x) 

H[cj>\ =£dzjd''-\ [i|V</)p + iro^2 ^ ^^^^4] (2.3) 

Parameter tq involves local random temperature fluctuations 5t(x) via to — mp + 5t(x), where {ST{x))conf = and 
{6T{x)6T{x'))conf = ^S{x — x') with A > 0. It should be noticed that ttiq corresponds to a "bare mass" representing 
linear temperature deviations from the mean - fleld critical temperature. 

In the present paper we wish to investigate the surface critical behavior of semi-inflnite Ising-like systems with 
quenched bulk randomness at the special transition. The critical behavior at the special transition has own pecu- 
liarities. In general case effective Hamiltonian for such systems must involve terms describing surface interactions 
@0j29lg|,|| 
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(2.4) 



Usually these external surface fields h are position and time dependent, so that means they may be used as sources. 
But in the present investigation we restrict ourselves to the case h — Q. Thus, the common form of the effective 
Hamiltonian describing the surface critical behavior of quenched dilute semi-infinite Ising-like systems in the replica 
limit n — > can be written in the form 

H{$) = dzj d''~\[\ I V0 p +\ml I $ p 

+ ^^oE0f + ^"0(10^] + ^ / d'-'rc,$\ (2.5) 
i=i ■ 

It should be mentioned that fields (/)(r, z) satisfy the Neumann boundary condition ||2^,^, so we have at z = that 
dz(l>{r, z) — 0. This Hamiltonian takes into account surface interaction in the form of additional term ^ J d'^^^rcocfP' ■ 
The model defined in ( ^.5[ ) is restricted to translations parallel to the boundaring surface. Thus, only parallel Fourier 
transformations in d — 1 dimensions take place. 



III. RENORMALIZATION OF THE CORRELATION FUNCTION 



The correlation function of the model (see Eq.(2_^)), which involves N fields 4>{'x.i) at distinct points Xi(l < i < N) 
in the bulk and M field 4){Yj, z — Q) = <f>s{Yj) at distinct surface points with parallel coordinates rj{l < j < M), has 
the form 



N M 

G(^'^'^)({xj{rj}) = (n0(xo n'^^(^^-)>- 



The corresponding full free propagator in the pz representation is given by 



G{p,z,z') 



1 
2ko 



Co + KO 



(3.1) 



(3.2) 



with the standard notation kq — -y/p^ + m^. Here, p is the value of parallel momentum associated with d — 1 
translationally invariant directions in the system. The first term in (3_^) corresponds to usual free bulk propagator 
in coordinate space, between the points r = (p, z) and r = (p, z ), and the second, so called 'surface' term, depends 
on the distance between the point r and its 'mirror image' r = (0, — z ). 

The formulation of the randomness problem in the spirit of approach introduced by Grinstein and Luther indicate 
that the renormalization process for the random systems is similar to the 'pure' case ]27| , |35t . As it is known p7| , ^ , 
in the theory of semi-infinite systems the bulk field 0(x) and the surface field (f>s{v) should be reparameterized by 
different uv-finite renormalization factors Zif,{u,v) and Zi{u,v) 



All 



and 



Mr) 



The renormalized correlation function involving N bulk and M surface fields {N, M) ^ (0, 2) can be written as 



(3.3) 



In order to obtain the critical exponent ryjj^ which characterizes surface correlations at special transition, it is 

sufficient to consider a two-point correlation function of surface fields G'^^''^\p) — {(fip, z ~ 0)if{—p, z =0)). 

It should be mentioned that the renormalized mass m, coupling constants u, v, and the renormalization factor 
Z^ are fixed via the standard normalization conditions of the infinite- volume theory |^,|^,^,^ . In order to remove 
short-distance singularities of the correlation function G^*^'^^ located in the vicinity of the surface and to define 
the surface-enhancement shift 6c and surface renormalization factor Zi, new normalization conditions should be 
introduced. We normalize the renormalized surface two-point correlation function in such a manner that at zero 
external momentum it should coincide with the lowest order of perturbation expansion of the surface susceptibility 
X||(p)-G(0'2)(p) 
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+ 0{uo,vo) 



Co + \/p^ + rri^ 

and its first derivatives with respect to . Thus we obtain necessary surface normahzation conditions 

G^'^^ (0; m, u, v, c) = — \- 
m + c 



and 



^G^'^'' {p; m^u,v, c) 



dp^ 



p=0 



2m{m + cy- 



(3.4) 



(3.5) 



(3.6) 



Eq.(3.5) defines the required surface-enhancement shift 6c. From Eq.(3.5) is easy to see that the surface susceptibihty 
diverge at to = c = 0. This point corresponds to tlie multicritical point (m^^^c^ ) at which special transition takes 
place. 

Ta king into account the above normalization condition (B.6) and expression for renormalized correlation function 
( |3.3|) it is possible to define the renormalization factor Zy = Z\Z^ via 



= lim^^[G("'2)(p)]-i. 
P^o p dp^ 



(3.7) 



It should be mentioned that the next singular behavior of these renormalization Z factors in the critical region 
takes place 



Z^ cx to'', 
oc to"! 



(3.8) 



where to is identified via the inverse bulk correlation length cx i'^, t — ^ rp " . Here ri is the standard bulk 

^ rl II li 

correlation exponent and exponent ry*^ is specific for our quenched random semi- infinite system. As it is known 35 
these exponents rj and rj^ arise as a RG arguments of an inhomogeneous Callan-Symanzik equation for correlation 

functions G^'^"* (3.3) 



^4' = rn-^ InZ^ \pp , 

^(0,2) 



d 

Vi' = "1— InZi I pp . 
dm 



(3.9) 



The simple scaling dimensional analysis of Gp" and mass dependence of Z factors (3.S) defines the surface correlation 
exponent rj^f via 



sp sp , 

V{ = Vi + V- 



(3.10) 



According to (3.7),(p.S|) and ( 3.10 ) we obtain for surface correlation exponent ryjj^ 

d 



m— lnZ|| 
om 



\FP 



91nZ||(w, u) d\iiZii{u,v) 
I3u{u,v) h f3v{u,v)- 



di 



dv 



(3.11) 



FP 



The above value should be calculated at the infrared-stable random fixed point (FP) of the underlying bulk theory. 
The other surface critical exponents of the special transition can be determined via the set of surface scaling relations 



IV. THE PERTURBATION SERIES UP TO TWO-LOOPS. 



As was indicated in the previous section, the calculation of the surface critical exponent t^JJ^ can be performed in 
accordance with Eq.( 3.11 ), where renormalization factor Z|| is defined via (3/7). It should be mentioned that here. 
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by analogy with the infinite- volume theory, we start from inverse surface correlation function [G^'^"*]"^ in order to 
avoid the external lines depending both on the external momentum p and the surface enhancement Cq in each external 
propagator 



G(p;z,0) 



Kq + Co 



(4.1) 



In other words such procedure of the correlation function inversion allows to amputate the denominators kq + cq 
from each external propagator. Thus we consider the Feynman diagram expansion of the unrenormalized surface 
correlation function [G'^*''^^]^^ in terms of the free propagator (3.2) up to the two-loop order 



[G^°-'^\p;mo,uo,vo)] ^ = Co + kq 



[cq + Ko\ 



e -0 



(4.2) 



Here the full internal lines denote the full free propagator (3.2) and the dashed external lines correspond to the factors 
g-fcoz_ accordance with the presence solely d—1 translational invariance in (2.5) only a partial cancellation in the 
external legs takes place and as a result the one-line reducible graphs do not disappear. All these graphs in (4.2) have 
their own weights arising from the standard symmetry factors of the effective Hamiltonian (2.5) 



with 



.(0) 



,(0) 



with t 



(0) 



3 

-2 



■Uq+Vq, 



■^0 + ^0+ 2woWo , 



AO) 



where — i J"-* corresponds to the one-loo p di agram, t^'' to the two-loop melon-like diagram, t'"^'' and t]^'' to the reducible 



and ^ with if =tf = (4"))^ 

,(0) 



(4.3a) 
(4.3b) 
(4.3c) 



(0) 



,(0) 



and irreducible two -loop diagrams in (O), respectively. With a view to avoid the usual bulk uv singularities, which 
are present in (4.2), we perform the mass reno rma lization (see [p6[). After the mass renormalization we obtain that 



the divergent parts of the correlation function (4^) associated with the bulk terms of the free propagator (3^) cancel 
and the "melon-like" diagram becomes subtracted. Thus we obtain the next expression for inverse correlation function 

[G(".2)i-i 



X 1 ^"--^ d 



2k ' 2k dk^ 
+ 



fc2=0 



G G 



(4.4) 



The superscript G means the full free propagator (3.2), its bulk and surface parts are denoted with "-" and 
signs, respectively. The line labeled by D represents the Dirichlet propagator 



G(p,z,z') 



1 
2ko 



which arises from the difference 



+ + 



Z... - [cq + Kq] 1 



(4.5) 



(4.6) 



The expansion (4.4) still holds the uv divergences which are connected with the presence of the surface in the 
system. These divergent parts represented by the closed lines with the index " +" . The short-distance singularities of 
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the inverse correlation function [G*-"'^^ (p)]^^ can be removed via the renormaliza tion of the surface enhancement. The 
surface enhancement renormalization is defined via the normalization condition ( |3.5| ). This condition can be rewritten 
for our inverse unrcnormalized correlation function in the form 



Z|| [G'^^'^^ (0; Too, Wo, ^^0, Co)] ^ = to + c. 



(4.7) 



Substituting [G(o^2)(p)]-i (ggg Eq.Q) at zero external momentum p = into the normalization condition (^J) gives 
us the equation for the surface enhancement shift 6c 



Sc = ^ — 1)(to + c) + cro(0; TO, Co = c + Sc). 



(4.8) 



This equation can be resolved using method of sequential iteratio n to determine the dependence of co = c + (5c on c 
and TO. Here fTo(0; cq) denotes the sum of all loop diagrams in (4.4). Among them cri corresponds to the one-loop 
graph, (72 denotes the melon-like two-loop diagrams 



(T2(p;to, Co) 



1 , 



2k " 



TO^ d 



(4.9) 



fc2=0 



(73 and (74 represent the last two terms in E(5j4.4), respectively. 

Taking into account Eq.(4.8) the expansion ( |4.4| ) For inverse surface correlation function can be written in the form 

(4.10) 



[G(°'^^(p;to, c)] ^ = K - TO + Z|| ^(to 4- c) - [cro(p;TO, co) - (7o(0; to, co)]. 



According to Eq.([4.10|), the surface enhancement renormalization of the inverse surface correlation function reduces 
to the subtraction from each diagram contribution the same diagram at zero external momentum. Thus, the uv 
divergences connected with the presence of the surface i n th e system are mutually canceled. 

The iteration process of finding solution of the Eq. (4.8) implies that at zero-order we put 5c — and find the 
first-order solution 



Sc^^'l = {Z^r^ - 1)(to + c)+ CTi(0; TO, Co = c). 



(4.11) 



Taking into account that we have dealt with the special transition where c — 0, for the renormalization factor Z|| 
according to Eq.(3/7) at one- loop order we obtain 



(Z|7i)(i) = - 2r 



5(71 (p; TO, 0) 



jio) 
'-I 



(4.12) 



p=0 



where ti = ti / {%-Km) . At two-loop order, the first-order contribution (5c^^^ must be taken into account in <Ji because 
in common we have 



(7i {p; m,CQ) = cTi (p; to, co) Iq 



9(71 (p; TO, Co) 



dco 



(5c(i) 



(4.13) 



The uv singular part of the second term of this expansion remove the divergence of the correspondent surface subgraph 
in (74, denoted via . As we can see from (4.4), the similar subgraphs arise twice in (73. But the contribution (73 is uv 
finite because the Dirichlet propagator Gd vanishes whenever z or z approach zero. Thus, the perturbation expansion 
of the inverse surface correlation function of our system up to two-loop order after bulk and surface enhancement 
renormalization can be written in the form 



[G(°^2)]-i (Z|71)(i)to - ((71 (p; TO, 0) - ai(0; to, 0)) -t- 



5(7i(0;to, Co) 



dco 



,5c(i) 



9(71 (p; TO, Co) 



dcn 



4 

^c(i) -^(7,(p;to,0). 

i=2 



(4.14) 



Following to the definition (B.7) we obtain the next expansion for the renormalization factor Zu up-to to loop-order 



Zu 



^?)^l 



1 - ^ + ^— (- - ln2) - lim 



TO d daiip; to, co) 



p^o p dp 



dcQ 



m d 

(7i(0;to, 0) - lim — — CTi(p; to, 0). 
I P^O p op ^—^ 



(4.15) 
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Performing the integration of Feynman integrals in ( 4.15 ) by analogy with [ p5[ we derive the result 



^0) 



4 2 / ^ ' 



where the constant A arise from the two-loop contribution in (4.15), 

A ~ 0.202428. 



(4.16) 



(4.17) 



Here the renormalization factor is expressed as a second-order series expansion in powers of bare dimensionless 
parameters Mq — Uo/(87rTO) and Vq = Vq/ (Sirm). The corresponding weighting factors iff^ and i^""* are obtained by re- 



placements (uo, vq) — > (mq, vq) in the original combinations t^^^ and fj'^'' from (4.3). As it is usual in sMperrenormalizable 
theories the next step is the vertex renormalizations 



.(0) 



n + 8 



6 



■u + V 



(4.18a) 



Again, the vertex renormalization at d = 3 is a finite reparameterization. The result for a modified series expansion 
of Z|| in terms of a new renormalized coupling constants u and w, normalized in a standard fashion (4.18a), is 



7-1.- -, n+2_ V 



1 1 



n+2 



A - i + ^^{ln^2 - ln2) ) 



12 



A - - + -{ln^2-ln2) v 



1 n+2 , 



2{ A - - + —^{ln^2 - ln2) ) uv. 



(4.19) 



Combining the renormalization factor Zy (u, v) together with the one-loop pieces of the beta functions 

n+8 _ 



Puiu,v) = -u(l 
I3v{u,v) = -vi^ 



6 



1 v-2u 

2 



through (3.11), we obtain the desired series expansion for rj^^' 



T?|P(u,u) = 

{n+2) 



+12- 



n+8)2""^'~' ■ 9 
where A{n) is a function of the replica number n, defined as 

A{n) ^2A + 



71+2 V 

"2(n+8) 6 

4 8 
A{n)u'^ + -A{l)v'^ H A{n)uv, 



(4.20a) 



(4.21) 



48 
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n+8 



{ln^2~ln2), 



(4.22) 



and renormalized coupling constants u and v, normalized in a standard fashion u=^^^^u and v=^v. 

In fact, the last expression (4.21) for i]^^^ provides a result for the cubic anisotropic model given by the effective 
Hamiltonian (2.5) with general number n of order-parameter components. In the case of infinite space, this cubic 
anisotropic model attracted much attention very recently (see e.g. [|l9|,^^ and references therein). 

In the present paper we restrict our discussion to the case of semi-infinite random Ising-like systems by taking the 
replica limit n ^ 0. According to (4.21) wc obtain the next two-loop expansion for the surface critical exponent 
at the special transition 



u w 3 2 4 , , 2 



+ -A{Q)u' + -A{l)v' + A{Q)uv. 
6 8 9 



(4.23) 



As it is well known, the knowledge of one surface critical exponent gets access via usual scaling relations to the 
other surface critical exponents. For convenience further below we suppress the superscript sp at the surface critical 
exponents. 
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V. CALCULATION OF THE SURFACE CRITICAL EXPONENTS 



The present section is devoted to numerical calculation of the surface critical exponents at the special transition. 
The individual RG series expansions for other surface critical exponents can be derived from (4.23) through standard 
scaling relations (with d = 3) 







2 ' 








= i^(l - v\\ 


), 




), 




), 


^(^-, 


A dH 


h2-77 


Pi d- 


2 + ' 


Ai 




/3i d 


+ 



V-L - 

711 

7i = K2-?7±), (5.1) 
Ai = 

<5i = 

Sn - 

Each of these critical exponents characterizes certain properties of the system with the surface in the vicinity of the 
critical point (see |^^). The values v, rj, and A = iy{d + 2 — ri)/2 are the standard bulk exponents. The correspondent 
series expansions for i/ and 77 at c? = 3 are given by ( ]^-p^ ) 



_ 1 
~ 2 

1 
324 
8 

~ 27 



(n + 2) 



11 . 
9 

27 



2{n + 8) 
2 



2uv 



(27n - 38)™ - ^^^^-^{27n - 38)w2 



[n + \ 



(n + 2) 2' 



3(n + 8) (n + 8) 



(5.2) 



For each of surface critical exponents we obtain from (5T) and ( 4.21 ) at d = 3 a double series expansion in powers 
of u and v truncated at the second order 



(5.3) 



],l>o 



Since perturbation expansions of this kind are generally divergent , the sufficiently powerful resummation procedure 
of the series is essential to obtain accurate estimates of the critical exponents. One of the simplest way is to calculate 
for each quantity we consider a sequence of rational Fade approximants in two variables from the original series 
expansions. This should work already well when the series behave in lowest orders 'in a convergent fashion'. Besides, 
if the series are alternating in sign |55[ , we can use more modern Pade-Borel resummation procedure | |5^ for their 
analysis. 

The results of our Fade and Fade-Borel analysis for the surface critical exponents at the special transition are 
represented in Table I. We evaluate the exponents at the standard RG random fixed point M 



-0.60509, 



2.39631. 



(5.4) 



The values [0/0], [1/0], and [2/0] are simply the direct partial sums up to the zeroth, first, and second orders, 
respectively. Fade approximants [0/1] and [0/2] represent the partial sums of the inverse series expansions up to the 
first and second order. 

Besides, by analogy with [ ^ we performed the calculation of nearly-diagonal two-variable rational approximants 
of the type 



[11/1] 



1 + aiu + aiv + aiiuv 
1 + biu + biv 



(5.5) 



and 
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1 + Oiu + biv + biiuv 

which are given in eighth and ninth columns of Table 1. 

The second column of Table 1 contains the ratios of magnitudes of first- (Oii) and second-order (02«) perturbative 
corrections appearing in inverse series expansions of our critical exponents. The larger (absolute) values of this ratios 
correspond to the better apparent convergence of truncated series. It is easy to see that the series of inverse expansions 
for all critical exponents, except (3i, are alternating in sign and consequently adapted to the above-mentioned Pade- 
Borel resummation analysis (see Appendix 1). Among the direct series the situation is more complicated. The ratios of 
first- Oi and second-order O2 perturbative corrections of the direct series expansions for the critical exponents 61, 711 
and 71 are positive [|60| . This means that the signs of the first- and second-order corrections do not alternate and hence 
the correspondent series are not suitable to the Pade-Borel resummation technique, since the [11/1] approximant of the 
Borel transform have a pole in the integration range. But these series are slowly convergent, because the contribution 
of the second-order are considerably less than contribution of the first-order. For example the ratio O1/O2 for the 
critical exponent Ai will be equal 35.1. Thus the above mentioned series adapted to the Fade analysis. It should be 
noted, that a very similar situation has been meet in the analysis of the surface critical exponents perturbation series 
expansions at the ordinary transition in pure js^ j and quenched dilute semi- infinite Ising-like systems | |3^ ]. 

The results of Fade-Borel analysis of the inverse series expansions are given in the last column of the Table I. 
These values give the most reliable numerical estimates. In common, we assume that the results obtained from the 
Fade-Borel analysis are the best we could achieve from the available knowledge about the series expansions in 
the frames of the present approximation scheme 

ry||= -0.238, Ai 1.101, ry^ = -0.116, /?i = 0.258, 

= 0.845, 71 = 1.442, Si = 6.343, Sn = 4.172. (5.7) 

Here the estimate of /3i, for which no Fade-Borel approximant R^^ exists, has been derived from the scaling relation 
= |((i - 2 + ?7||), where ^ = 0.678 [| and the above value of rjn = -0.238. 

The deviations of estimates ( |3.7| ) from the other second-order estimates of the table might serve as a rough measure 
of the achieved numerical accuracy. 



VI. EXPANSION 

As was mentioned above, there is an alternative method to analyze the influence of randomness on the critical 
behavior based on the renormalization-group approach introduced by Harris and Lubensky [Q. This method was 
employed by Ohno and Okabe p5[ for study critical behavior of semi-infinite systems with a Gaussian randomness 
in 4 — e dimensions in the frames of ^/e expansion for obtaining the two-loop approximation for correlation function 
and deriving corresponding series expansions for the surface critical exponents ?7|| and ri±. Their results in the case 
of special transition at n = 1 (see |^5|) with corresponding changes of coupling constants normalizations {u — > u/24, 
w — m/3) in accordance with our notations are written in the form 

u V 7 o 11 o 7 ^,3.. 

7711 = \ u -\ V + - uv + 0(e2) 

'II 3 2 12 12 4 ^ ^ ' 

u w 11 9 23 , 11 ^,3. ,„ ^. 

.. = -g-44--u^ + -.^ + -.. + 0(..). (6.1) 

Unfortunately, the surface critical exponents r/|| and r/x have been obtained only to the first order in ^/e from these 
equations in |^5| . In the present paper we derive the next term in -^e expansion for above mentioned surface critical 
exponents using the fixed-point values up to 0(e) |^,|^ 

* „ , , ^^ 110 + 63C(3) 
" =-'^^53+^' 53^ ^' 




where C(3) ~ 1.2020569 is the Ricmann ^-function, and the usual geometric factor Kd = /T{d/2) has been 

absorbed into the redefinitions of the coupling constants. As a result we obtain 
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^11 




756C(3) - 641^ 



2 • 53 



2 



378C(3) - 347 
^^^ = 'Vl06 + 2-53^ " 

Taking into account scaling relations for surface critical exponents and i/e expansions for random bulk exponents 
V and r\ we obtain respective perturbation series expansions for other surface critical exponents 




16-532 
3 (252C(3) - 461) 
16 532 
2268C(3) - 2453 
8^532 ' 
3 (378C(3) - 347) 
532 ' 



2268C(3) - 

711 = o + 7\/^ 

71 = 



, (630C(3)- 1073) 
5i = 5 + 5,/--3 

. ^ (756C(3)- 1277) 
- ^ + V 53 - 2 53^ " (^-'^ 

Similarly as in the case d = 3 of the previous section, we perform a Pade analysis of our expansions at e = 1 . The 
numerical values of surface critical exponents obtained in this way are represented in Table II. It should be noticed 
that as was shown in [[40y4l| ], the \fi expansion is not Borel summable. 

The Pade approximants [1/0] for the exponents 7711 and r/^ reproduce the first-order results obtained by Ohno and 
Okabe [Q. On the other hand, the other exponents, /3i, 711 and 71 slightly differ from the previous results |^ 

/3i = 0.17, 711 = 0.78, 71 = 1.26. (6.5) 

The reason is that we calculated our [1/0] estimates directly from each respective -^e expansion, while in Ref. [|5| they 
were obtained from the scaling relations using the above mentioned numerical vales of r^n and 77^ ( |6.3| ). In addition we 
performed the analysis of the correspondent series expansions for the surface magnetic shift exponent Ai, exponents 
8\ and 8\\ which give relations between the surface magnetization and the surface and bulk external magnetic fields, 
respectively. 

It can be easily verified that the above first-order approximants of the critical exponents satisfy the scaling relation 

= ('/ + ^||)/2 (6.6) 

with the value of the bulk exponent r\ — — + O(e^) 

Comparing the results given in Tables I and II we see that the values of the first-order approximants, denoted by 
[1/0] and [0/1] in both cases, are of comparable magnitudes. But, on the other hand, the values of the second-order 
approximants are significantly different in both tables. As was shown previously [|37|-|4l[], -^e series expansions possess 
rather irregular structure and are practically unsuitable for subsequent resummation and are ineffective for obtaining 
relia ble n umerical estimates. Our results confirm this assumption. If we try to reproduce the numerical value of 77 
(see ( |6.6| )) from our second-order data of Table II, we always obtain negative values. But, this does not agree with the 
sufficiently precise positive results of massive field-theoretic approach for random bulk systems in three dimensions up 
to two-loop H,^, three-loop [ p7|pO|j40[ ] and to the four-loop [|5^ order. Besides, very recently was obtained the value 
f] = 0.025 ±0.01 in the frames of five-loop renormalization-group expansions [|l^ and up-to six-loop order 77 — 0.030(3) 
[Q). This discr epan cy is not present in our calculation performed directly at d = 3 (see Table I). From the surface 
scaling relation (|6.6D and the second-order results of Table I we obtain the value 77 = 0.031, which quite well agree 
with previous estimates. 

VII. SUMMARY 

The main aim of the present article was an investigation of the influence of quenched bulk randomness on the 
surface critical behavior of semi-infinite Ising-like systems at the special transition. To summarize, we have calculated 
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the surface critical exponents of the special transition for such systems using two alternative technique: the massive 
field qtheory approach directly at d = 3 dimensions up-to two-loop order, and the ^fl expansion at d = 4 — e 
dimensions to the order of 0((-ye)^). In the last case we extend up to the next to leading order, the previous first- 
order results obtained by Ohno and Okabe |Q. But, the -^/e expansions possess rather irregular structure, as was 
shown in This makes them practically unsuitable for subsequent resummation and ineffective for getting 

quantitative numerical estimates. 

In both cases the resummation of obtained perturbation series expansions for surface critical exponents was per- 
formed using Pade analysis. Moreover, in the case of the massive field theory the resulting two- loop series expansions 
were resummed by means of more precise Pade-Borel resummation technique. In the previous section we have dis- 
cussed some merits of using the massive field-theoretic approach directly in d = 3 dimensions for obtaining most 
accurate numerical estimates for critical exponents. Thus, the best estimates for surface critical exponents of semi- 
infinite systems with quenched bulk disorder at the special transition, which we can obtain in the frames of the present 
approximation scheme, are 

?7||= -0.238, Ai = 1.101, = -0.116, /3i = 0.258, 

711 = 0.845, 71 = 1.442, (5i = 6.343, hi = 4.172. (7.1) 
The obtained results are evidently different from results obtained for pure semi-infinite Ising-like systems [^ , |6l| 

?7|l = -0.165, Ai = 0.997, 7/^ = -0.067, /3i = 0.263, 

711 = 0.734, 71 = 1.302, (5i = 5.951, = 3.791 (7.2) 

and confirm the assumption that the presence of quenched bulk disorder affects the critical behavior of boundaring 
surface. So, in the case of special transition, similarly as in the case of previously investigated ordinary transition 
p6| , the new set of surface critical exponents appear. It should be mentioned that at the present time the value of 
the surface crossover exponent 4> for semi-infinite systems with quenched bulk disorder is still an open question. This 
question will be the topic of our next publication. 
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APPENDIX 



As was mentioned above, the power series expansions of surface critical exponents (see Eq.( |5.3D ) are in general 
not convergent. In order to obtain meaningful and rather accurate numerical estimates we must apply to them 
sufficiently powerful "resummation" procedure. In the present paper we employ a two- variable resummation technique 
P, ^7p^j5^ , |59| | which is a simple generalization of the sing le- variable Pade-Borel method. The starting point of this 
calculation is construction for truncated power series (^.3[) the Borel transform 



F{utM) - 51 



+ 



Then we construct the rational approximant F^{x, y) 



F^{u,v) = 



1 + aiQU + aoiv + aiiuv 
1 + biou + boiv 



(7.3) 



(7.4) 



which is extrapolation of the Borel transform (L3). It is clear that at w = or t; = we obtain from (7^) the usual 
[1/1] Pade approximan t. T he coefficients aji and bji are expressed via the appropriate expansion coefficients of the 
initial function f{u,v) (5.3) 



o-oi — .910 + ^10, ^10 = —920/ 910, 
o-oi = .901 + ^01, &01 = —302/501 J 
an = 911 + ^10501 + ^oiffio, 



(7.5) 
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where gji = fji/{j + 1)1. Hence, for the resummed function by means of the Pade-Borel resummation technique we 
obtain 

f{u,v)= F^{ut,vt)e-^dt. (7.6) 
Jo 

This technique has been apphed to the direct power series expansions of surface critical exponents to find the corre- 
sponding Pade-Borel approximants R~^. 

TABLE I. Surface critical exponents of the special transition for d = 3 up to two-loop order at the random-fixed point 
u* = -0.60509, V* = 2.39631. 



exp 


On 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[11/1] 


[1/11] 




m 


-23.82 


0.00 


-0.324 


-0.245 


-0.205 


-0.237 


-0.244 


-0.238 


-0.238 


Ai 


-3.39 


0.75 


1.074 


1.229 


1.083 


1.046 


1.083 


1.090 


1.101 


V± 


-3.35 


0.00 


-0.162 


-0.139 


-0.087 


-0.102 


-0.115 


-0.114 


-0.116 


Pi 


0.00 


0.25 


0.25 


0.25 


0.263 


0.263 








711 


-3.14 


0.50 


0.824 


0.979 


0.825 


0.783 


0.825 


0.834 


0.845 


71 


-2.56 


1.00 


1.405 


1.680 


1.410 


1.327 


1.410 


1.421 


1.442 


5i 


-1.41 


5.00 


6.619 


7.394 


7.062 


5.521 


6.205 


6.236 


6.343 


Sii 


-1.40 


3.00 


4.295 


5.279 


3.926 


3.418 


4.032 


4.070 


4.172 



TABLE IL Surface critical exponents of the special transition from the ^/e expansion. 



exp 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 


m 


0.00 


-0.3.36 


-0.252 


-0.289 


-0.287 


-0.295 


Ai 


0.75 


0.960 


1.016 


0.938 


0.917 


0.940 


V-L 


0.00 


-0.168 


-0.144 


-0.149 


-0.151 


-0.151 


Pi 


0.25 


0.208 


0.210 


0.197 


0.198 


0.194 


711 


0.50 


0.752 


0.838 


0.740 


0.714 


0.741 


71 


1.00 


1.252 


1.338 


1.224 


1.190 


1.227 


5i 


5.00 


6.682 


7.535 


7.019 


6.729 


7.104 


Su 


3.00 


4.346 


5.441 


4.608 


4.232 


4.671 
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